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Abstract. The ground state correlations induced by a general pairing Hamiltonian in 
a finite system of like fcrmions are described in terms of four-body correlated structures 
(quartets). These are real superpositions of products of two pairs of particles in time- 
reversed states. Quartets are determined variationally through an iterative sequence 
of diagonalizations of the Hamiltonian in restricted model spaces and are, in principle, 
all distinct from one another. The ground state is represented as a product of quartets 
to which, depending on the number of particles (supposed to be even, in any case), 
an extra collective pair is added. The extra pair is also determined variationally. 
In case of pairing in a spherically symmetric mean field, both the quartets and the 
extra pair (if any) are characterized by a total angular momentum J = 0. Realistic 
applications of the quartet formalism are carried out for the Sn isotopes with the 
valence neutrons in the 50-82 neutron shell. Exact ground state correlation energies, 
occupation numbers and pair transfer matrix elements are reproduced to a very high 
degree of precision. The formalism also lends itself to a straightforward and accurate 
description of the lowest seniority and 2 excited states of the pairing Hamiltonian. 
A simplified representation of the ground state as a product of identical quartets is 
eventually discussed and found to improve considerably upon the more traditional 
particle-number projectcd-BCS approach. 
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1. Introduction 

Pairing strongly affects the behavior of a large variety of quantum many-body systems 
ranging from very small, like nuclei, to very large, like stars [lj. The concept of pairing 
made its first appearance in the condensed matter context in an attempt to explain the 
phenomenon of electron superconductivity in metals [2 [3]. Soon after, this concept was 
imported in nuclear physics [I] and, since then, it has proved crucial to understand a 
whole series of major properties in nuclei such as the odd-even staggering in binding 
energies and the presence of an energy gap in the spectrum of even-even nuclei [5] [6] . 
In recent years, in conjunction with an increased experimental exploration of nuclei far 
from stability, a revival of interest in pairing has been observed due to its fundamental 
support to the theoretical description of loosely bound nuclear systems [IJ |5]. 

In a mean-field treatment of nuclei, pairing mimics the short range part of the 
effective nucleon-nucleon interaction. Its basic action on nucleons is that of scattering 
them between pairs of time-conjugate single-particle levels. In its simplest form, 
namely assuming this scattering to be independent of the levels involved, the pairing 
Hamiltonian reads as 

H = Y,tkN k -gY / PlPk>, (1) 

k kk' 

where Pi creates a pair of particles in time-conjugate states on the level k of the mean 
field, N]~ counts the number of particles on this level, is the associate energy and g 
the constant pairing strength. In spite of being highly schematic, this Hamiltonian has 
been applied to various fields of physics ranging from nuclear structure to condensed 
matter, where its applications include not only macroscopic superconducting systems 
but also finite size microscopic systems such as ultrasmall metallic grains [7]. 

The Hamiltonian (TjQ) is commonly treated in the BCS approximation [3] or in 
the more refined Hartree-Fock-Bogoliubov (HFB) approach |6J. As it is well known, 
these methods, while suited to describe macroscopic systems, being exact in the 
thermodynamical limit [12], show some limitations when applied to systems with a 
limited number of particles such as nuclei. This is mainly due to fact that they make 
use of a wave function that is not an eigenstate of the number operator. 

In a series of papers in the 60's, Richardson [HI El HD], also in collaboration with 
Sherman studied the pairing Hamiltonian (TfJ) and showed that the corresponding 
eigenvalue problem can be solved exactly in a semi-analytical way. More precisely, 
it was shown that, for an even number of particles (2iV), the exact eigenstates of the 
Hamiltonian (1) can be written as a product of distinct, collective pairs of particles built 
on time-conjugate states. These pairs can be either real or complex, with reference to 
their mixing amplitudes, depending on the value of the strength g. Complex pairs, 
in particular, always occur in a complex-conjugate form. In [13], it was shown that 
a product of two such complex-conjugate pairs actually forms a four-body correlated 
structure with real mixing amplitudes (quartet). The exact ground state of the 
Hamiltonian (CQ) can therefore be formulated in an equivalent form as a product of 
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these quartets. Quartets are expected to be all distinct from one another and, for 
odd values of N, they are joined in by an extra collective pair also characterized by 
real amplitudes. Such an exact and (semi) analytic representation of the ground state 
in terms of quartets holds, however, only for the constant pairing Hamiltonian (1). 
Extending this formalism of quartets to a general pairing Hamiltonian and investigating 
how much this new representation improves upon more traditional descriptions of pairing 
are the basic motivations of the work presented here. 

Quartets will be constructed via an iterative variational approach based on 
diagonalizations of the Hamiltonian in restricted model spaces. A similar technique, but 
based on pairs, has been employed in a previous work on pairing [H] and a comparison 
with its results will be discussed below. Each diagonalization is meant to update a 
quartet while driving the ground state towards its minimum. Whenever needed, an 
extra real collective pair, also determined variationally, will be added to these quartets. 
As a testing ground for the present approach, we will examine the Sn isotopes with the 
valence neutrons in the 50-82 neutron shell. Comparisons will be carried out with the 
exact solutions as well as with more traditional approaches to pairing. The analysis will 
not be limited to the ground state only but it will concern low-lying excited states as 
well. As a final step, a simplified representation of the ground state as a product of 
identical quartets will be formulated and compared with its direct pair analogue, i.e. 
the particle-number projected-BCS ground state [6J. 

The paper is organized as follows. In section 2, we review the Richardson approach 
to the constant pairing Hamiltonian and point out the emergence of quartets in its 
ground state. In section 3, we describe the procedure which has been developed to 
construct the quartets (and the extra pair) in the case of a general pairing Hamiltonian. 
In section 4, we present the results relative to the ground state while, in section 
5, we extend the quartet formalism to excited states. In section 6, we discuss the 
representation of the ground state as a product of identical quartets and, finally, in 
section 7, we summarize the results and draw the conclusions. 



2. Quartets in the Richardson formalism 

This section will be devoted to an analysis of the constant pairing Hamiltonian 

H = Y^^-gJ2 L l L ^ ( 2 ) 

i ii' 

where 

•Ni — a imj a imj ; = -fjimj) ^irrii = a imfi , im i - (3) 

mi mi>0 

The operator a\ m . creates a particle on the single particle level i, with energy 
and degeneracy a\ m . being the time reversed operator. In the case of pairing in 
a spherically symmetric mean field, i = Ui,li,ji (adopting the standard notation), 
Qi = 2ji + 1 and a\ m% = (— l) Ji ~ m *af_ , raj being the projection of In this case, 
the operator L\ creates a pair of particles on the level i with total angular momentum 
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J = 0. In the absence of spherical symmetry, L\ simply creates a pair of particles with 
opposite projections of the intrinsic spin on the doubly degenerate level i (m 8 = ± and 
Qi = 2 in this case). 

By further defining 

Ni mi = dim^imi ~\~ di—mfii—rrni (4) 

the Hamiltonian (j2J) reduces to 

h= E tiNimt-g E E pLa^ (5) 

i,mi>0 «,mj>0 i',m i />0 

which, assuming the simplified notation = i,m, > 0, is nothing but the Hamiltonian 
([T]). The index /c will be herewith assumed to vary in the interval (1, Q). 

The eigenvalue problem relative to the Hamiltonian fl2]) can be solved by applying 
the Richardson formalism [8j [9j QUI QT]. One finds that, for a system of N pairs and by 
neglecting the presence of unpaired particles in the model space (these do not partecipate 
in the pair-scattering generated by the Hamiltonian (CQ) (TTJ), any eigenstate can be 
formulated as the pair product state 

n n i 

l*> = II s Jl°>> B i = y Lir=r p l ^ 

v=\ k=l Z€k Cl/ 

if the N parameters e v (the so-called pair energies) are roots of the set of iV coupled 
non-linear equations 

n a N 2a 

i-E— — + y -^—=o. (7) 

^ Of , — p ^ p — e 

The eigenvalue associated with \^) is just the sum of the corresponding pair 

energies, i.e. 

N 

£ W = E e - (8) 

Solving numerically the equations © meets some technical difficulties. Indeed, by 
varying the strength g, it can happen that two (real) pair energies e u become equal 
thus giving rise to a singularity. When this occurs, these energies turn from real into 
complex-conjugate pairs. It is in order to handle only real quantities that Richardson 
introduced some new variables £a and rjx [TTJ. In the case of the ground state and for 
an even number iV of pairs, these are defined such that 

e 2 A-i = 6-^A, e 2X = £ x + ir, x (A = 1, 2, ■ • ■, N/2). (9) 

£a is assumed to be always real while r]\ can be either pure imaginary (corresponding 
to real and distinct e2A-i and e2\) or real (corresponding to complex-conjugate pair 
energies). Once the transformation (jUJ) is made, the equations ([7j) depend only on the 
real quantities £\ and r]\ which can therefore be calculated numerically (although at the 
cost of a further change of variables [H]). With increasing g, the quantities r}\ turn one 
by one from negative to positive values so generating the complex-conjugate pairs. The 
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case of odd N can be described in the same formalism of equation Q by assuming that 
one pair energy is real and the remaining N — 1 pair energies occur in complex-conjugate 
pairs [TT] . 

In the new variables (9) the ground state for a system of an even number of pairs 
can be written as a product of complex-conjugate collective pairs 

N/2 

\* 9 s)= II^A-I^aIO). (10) 
A=l 

The complex-conjugate form (|HD guarantees that the product B^^B^ can be set into 
a real form no matter what g. Indeed, one finds that [13] 

^ A _i4A = (ri) 2 + ^(©i) 2 , (ii) 

where the pairs 

n = E r2£ 2 %"i A + „ 2 f^ETiM, (12) 

fc=i l^ e fc - f a) + ?7a fc=i 

have amplitudes 7^ and 0^ depending only on the real coefficients £a and r/^. More 
explicitly, one has 

SiL-xBl = E (T^# + rj&tHW)PtPl (14) 
fc,fc'=i 

This expression holds true both for real and complex-conjugate pairs. The product 
Blx-iBlx can therefore be always represented as a real superposition of products of two 
pairs of particles in time-reversed states. In the following, we will refer to a four-body 
correlated structure having these general properties as a quartet. The ground state of 
the pairing Hamiltonian for a system with an even number N of pairs consequently 
acquires the form of a product of quartets: these are, in general, partly amenable to 
the product of two real collective pairs and partly associated with the product of two 
complex-conjugate pairs. For a system with odd N, instead, the ground state can always 
be reduced to a product of quartets and an extra real collective pair. 



3. Quartets for a general pairing Hamiltonian 

The analysis of section 2 cannot be extended to the more general pairing Hamiltonian 
H = Y J £iMi-Y J gwL\L i , (15) 

i ii 1 

since no analytical solution of the eigenvalue problem is known in this case. It is 
nevertheless possible to carry out an approximate treatment of the eigenstates of 
Hamiltonian (15) in terms of quartets by resorting to numerical methods. In this section, 
we will show how this has been realized. 
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We begin by considering a system with an even number N of pairs and we assume 
a ground state of the form 

n/2 (i,n) 

l*p-> = II Qt\o), Ql=H QwHpI- (16) 

u=l k<k' 

The operator Q* v is a generalization of the operator (|14p and creates a quartet with real 
amplitudes q^, to be determined. \^ gs ) is therefore a product of quartets which are, in 
principle, all distinct from one another. In order to search for the most appropriate q^i's 
we make use of an iterative variational procedure that resorts to diagonalizations of the 
Hamiltonian in spaces of a rather limited size. This procedure draws inspiration from 
an analogous technique previously developed for a treatment of pairing correlations in 
terms of a set of independent pairs [H] and works as follows. Let us suppose that, at 
a given stage of the iterative process, one knows the state (fTB|) and let us construct the 
space 

( N/2 \ 

F {p) = \plPl n Qlm ■ (17) 

I u =l(v^p) ) l<k<k'<n 

The states of F^ are generated by acting with the operators P^Py (k < k') on the 
product of all the quartets Q\ but the p-th one. The dimension of each space (TTTj) is 
therefore Q(Q — l)/2 and one can form at most N/2 such spaces. By diagonalizing the 
Hamiltonian in and searching for the lowest eigenstate, one constructs the state 

N/2 

l¥ ne w)) = Q Kne W ) J| Q t) >. (18) 

This differs from \^ gs ) only for the new quartet Q^ new ^ and its energy is by construction 
lower than (or, at worst, equal to) that of \^ gs ). As a result of this operation, then, the 
quartet Qji new ) has updated Qfi while all the other quartets have remained unchanged. 
At the same time the energy of \^ gs ) has been driven towards its minimum. Performing 
a series of diagonalizations of H in F^ for all possible p values (1 < p < N/2) exhausts 
what we define an iterative cycle. At the end of a cycle all the quartets Ql have been 
updated and a new cycle can then start. The sequence of iterative cycles goes on 
until the difference between the ground state energy at the end of two successive cycles 
becomes vanishingly small. 

The case of an odd number iV of pairs proceeds along the same path but one now 
assumes that the ground state is formed by (N — l)/2 quartets and one extra pair. One 
has, then, 

(n-i)/2 n 

hu=£ f n Gti°>. B^ = j2P k pi (19) 

v=l k=l 

the quartets Q\ being defined as in equation ([TBI) and the amplitudes (3 k being supposed 
to be real. In this case, the iterative cycle consists of two different steps. In the first 
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one, one performs a series of diagonalizations in the spaces 




v=i(v^p) ) i<k<k'<n 



(20) 



for all possible p values. These diagonalizations update the (N — l)/2 quartets while 
leaving unchanged the pair. Once these have been completed, as a second step, one 
performs a diagonalization in the space 



which updates the pair with the quartets acting as spectators. Also in this case, the 
sequence of iterative cycles goes on up to the convergence of the ground state energy. 

4. Ground state results 

As an application of the procedure just described, we have considered a case of pairing 
in a spherically symmetric mean field. This application has concerned, in particular, 
the Sn isotopes with the valence neutrons in the 50-82 neutron shell. As customary 
for these nuclei, the model space has been restricted to the five neutron orbitals 57/2, 
^5/2, ^3/2, S1/2, and /in/2. Single-particle energies ej and pairing strengths have 
been taken from the work of Zelevinsky and Volya [15] and are the same used in a 
previous analysis of pairing |14j . Within this model space, exact calculations have been 
performed by resorting to the method of [16] which is based on a classification of the 
many-body configurations within the seniority scheme p21 (THJ [19] and on the use of the 
SU(2) quasispin algebra. 

The quartets Q\ of equation (fTB"l) are defined in an uncoupled formalism and so 
they are not characterized by a good angular momentum J. The same is true for the 
pair of equation (fl9"j) . In spite of that, J = pair and quartets can be generated 
quite naturally by proceeding step by step for increasing values of N. For N — 1, the 
diagonalization of the Hamiltonian in the space S (j2"Tj) generates a J = pair as ground 
state. For iV = 2, a similar diagonalization in gives rise to a J = quartet. Moving 
to iV = 3 (a system which, in our formalism, is described as a product of one pair and 
one quartet), it is sufficient to use the previously generated J = pair and quartet as 
an initial ansatz to generate, at each step of the iterative cycle, either a J = pair or a 
J = quartet. The same mechanism holds true going to larger iV's. We notice that, as 
a result of the degeneracies of the Hamiltonian, the coefficients q^, that are generated 
are independent of the angular momentum projections and so the quartets (fTBl) acquire 
the simple form 



with the operators L\ <^ creating a J = pair on the orbital %. These J = quartets 
are therefore a linear combination of products of two uncorrelated J = pairs. 




(21) 




(22) 
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In figure 1, we have plotted the relative error in the ground state correlation 
energy (i.e., the energy of the ground state relative to the lowest uncorrelated many- 
body configuration) as a function of the number N of pairs. Results obtained with 
the procedure of section 3 are shown at various levels of approximation, each level 
corresponding to a different number of iterative cycles (indicated by the number next 
to each line). Already at the lowest level of approximation (1 cycle), the approach is 
seen to provide a remarkably good agreement with the exact results. This agreement 
progressively improves with increasing the number of iterative cycles up to becoming 
basically exact. 

It is of interest to compare these results with those which can be obtained in a 
formalism of independent J = pairs. The dot-dashed line in figure 1 refers to a 
calculation in which the pairing ground state has been assumed to be a product of real, 
collective, distinct, J = pairs [13] . The approach based on quartets appears to be 
considerably more effective than the corresponding one with J = pairs. 

In figure 2, we have plotted the root mean square values of the re lative errors in 
the occupation numbers, i.e. the quantity cr(oc) = \J (J2i=i (n))/5 with Ai(n) = 
^(exact) _ n (appr)y n (exact) ^ ^ e i nc i ex { running over the 5 orbitals of the shell 50-82 and 
rii being the corresponding occupation numbers. In figure 3, on e finds the ana logous 
rms values for the pair transfer matrix elements, i.e. a(tr) = J (X)f =1 Af (f))/5 where 

A,(t) = (tt act) - tt VVT) )lt { r Ct) and U = \(S>(N)\lI\V(N ~ The notation is the 
same as in figure 1, i.e., the dot-dashed line in both figures refers to the calculation 
with J = pairs while the remaining lines show the results obtained in the formalism 
with quartets at various levels of approximation. The results of these figures confirm 
the analysis of figure 1. 



5. Excited states 



The formalism of quartets illustrated in section 3 for the ground state also lends itself 
to a straightforward extension to excited states. To such a purpose, once the procedure 
relative to the ground state has been completed, we proceed by diagonalizing the 
Hamiltonian either in the space 

f N/2 \ 

G {E) = \PlPl II <#|0> (23) 

I v=l(u^p) ) X<k<k'<Xl,p=l,2,-,N/2 

or in the space 

t (JV-l)/2 (N-l)/2 \ 

G {o) = \pMrf n qvp)A n aim > ^ 

{ v =\{y±p) v=\ ) l<k<k'<n,p=l,2,-,(N-l)/2 

depending on the number of pairs iV: G^ E ' is employed for systems with even iV (i.e., 
in the presence of quartets only), when iV is odd. G^ includes at once all the 
states of the spaces FW (HZJ). These are, more specifically, all the states that can be 
formed by replacing, one at a time, each quartet of the ground state with all its basic 
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components P^fe'- * s formed, instead, by all the states of the spaces 

plus those of S (T2TT) . In addition to those states which are obtained by replacing each 
quartet with its components PlPy, therefore includes the states which result from 
the replacement of the pair with its components P\. One should notice that, as a 
result of the diagonalization of the Hamiltonian in these spaces, the excited states are 
no longer a simple product of quartets (and a pair, in case) as for the ground state but 
they are rather a linear combination of these products. These diagonalizations also lead 
to a rederivation of the ground state: this fact guarantees the orthogonality between 
this state and the excited ones. 

In tables 1 and 2, we show some examples of excitation spectra that are generated 
with this procedure. We adopt the same formalism used in the derivation of the exact 
eigenstates (section 4) and classify the approximate excited states according to the 
seniority quantum number v [I7J HEJ [19] . that is to say the number of nucleons not 
pairwise coupled to angular momentum zero. The energies shown in table 1 are those 
of the lowest 10 v — excited states in two middle shell Sn isotopes, namely 116 Sn and 
118 Sn. All these states are therefore characterized by a J = angular momentum. Table 
2 shows, instead, the energies of the lowest 10 v — 2 states still for 116 Sn and 118 Sn. In 
this case, each energy is usually associated with different J ^ degenerate eigenstates. 
The agreement between exact and approximate results is pretty good. 

Some comments on the internal structure of the excited states just discussed are in 
order. Quartets forming v = eigenstates are always of the type ff22l (J = quartets). 
In v = 2 eigenstates, instead, each component is characterized by one J ^ quartet of 
the type 

[«}«}]o(£^I)> (25) 

i 

the remaining ones being still J = quartets. This quartet is a product of a collective 
J = pair and an uncorrelated J ^ pair. As it can be seen from the above expression, 
following the general form ( fl6l) of the quartets, the latter pair is characterized by two 
nucleons with the same angular momentum j. We finally remark that these quartets are 
such that one cannot form v = 4 (or higher) states with good total angular momentum. 
The present formalism is therefore not well suited to describe these states. 

6. The ground state as a product of identical quartets 

The description of the ground state that has been pursued so far is based on quartets 
which are all distinct from one another. Treating systems of like particles (neutrons 
in the case just analyzed of Sn isotopes), quartets have been defined in terms of 
these particles only. In an attempt to describe proton-neutron pairing, use has often 
been made in literature of a-like quartets, namely four-body structures composed of 
two neutrons and two protons. In one of the earliest works on this subject [20J, 
a-like quartets have been introduced in a BCS-type wavefunction. Later on, more 
sophisticated wavefunctions have been proposed, still of BCS-type, which include the 
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contribution of quartet correlations to the proton- neutron pairing (e.g., [ZD [22]). As we ^ 
known, however, the use of BCS-type wavefunctions has the drawback of not conserving 
neither the particle number nor the isospin of the nucleus. Quite recently, an accurate 
description of the proton-neutron isovector pairing in nuclei has been obtained with a 
quartet model which does not violate any symmetry of the Hamiltonian j23j[24J. In this 
approach, the nuclear ground state is simply described as a product of identical a-like 
quartets with angular momentum J = and isospin T = 0. The present analysis offers 
the interesting opportunity to investigate, in a case of pairing among like particles, how 
well a description of the ground state as a product of identical quartets can approximate 
a description in which quartets are left free to be distinct from one another. To such 
a purpose, we have therefore compared the results of section 4 with those which are 
obtained by representing the ground state as 



in cases of odd N. In these expressions, we have kept the same definitions of Ql and B> 
given in equations (fT6l) and (fT9l) . respectively. The amplitudes q£y and (3k have been 
determined by minimizing the ground state energy. In figure 4, we compare the relative 
errors in the ground state correlation energy obtained in the approaches with identical 
quartets (triangles) and with independent quartets (squares). For completeness, in the 
same figure (circles), we also show the same quantity calculated by approximating the 
ground state as a product of identical pairs 5' only, that is to say in the well known 
particle-number projected-BCS approximation (PBCS) [6]. Clear differences can be 
observed between these three approaches, the corresponding errors often differing from 
one another by something like two orders of magnitude. The worst approximation 
appears to be the PBCS one (it should be noticed, however, that even in this case the 
relative errors remain confined within a quite acceptable range). The approach with 
identical quartets performs definitely better than the corresponding one with pairs and 
shows itself as a very good approximation of the exact ground state. 

7. Summary and conclusions 

In this work, we have provided a description of the low-lying spectrum of a general 
pairing Hamiltonian for a finite system of like particles in a formalism of quartets. These 
are, in general terms, real superpositions of products of two pairs in time-reversed states. 
Quartets have been first formulated in an analytical form in the case of a constant pairing 
Hamiltonian. This concept has been then extended to a general pairing Hamiltonian. To 
such a purpose, a numerical procedure has been proposed which is based on an iterative 
sequence of diagonalizations of the Hamiltonian. This procedure has been applied to 
the Sn isotopes with the valence neutrons in the 50-82 neutron shell. 




(26) 



in cases of even N, or as 




(27) 
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Two major features of this approach are worthy being emphasized: it does not 
violate any symmetry of the Hamiltonian and it reduces the original shell model 
diagonalization, usually involving very large spaces (up to 6x10 s basis states, in the 
uncoupled formalism, in the case of Sn isotopes [16]), to a sequence of diagonalizations 
in reduced spaces of only a few hundreds of states. As a major result, we have found 
that a product of distinct J = quartets (and, in case, an extra J = pair) can provide 
a basically exact representation of the pairing ground state. The quartet formalism 
has also been extended to excited states showing that it describes well the seniority 
and 2 eigenstates. Comparisons with a traditional approach to pairing like the particle- 
number projected-BCS approach, in which the ground state is supposed to be a product 
of identical J = pairs, or with a more sophisticated approximation based on distinct 
J = pairs, have evidenced a much higher quality of the approximation based on 
quartets. Finally, we have shown that a simplified description of the ground state as a 
product of identical quartets provides an excellent approximation of the exact ground 
state. 

The formalism of quartets developed in this work has shown to be quite successful 
in describing like-particle pairing in a realistic context. In order to be extended to the 
treatment of more general systems, such as protons and neutrons interacting through 
isovector and isoscalar pairing forces, this formalism needs to be further developed by a 
redefinition of the basic operators forming the quartets. This development is the scope 
of future studies. 
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Table 1. Comparison between exact and approximate energies of the lowest 10 v = 
excited states in 116 Sn and 118 Sn. Approximate results have been obtained with the 
procedure of section 5. All values are in MeV. 



116g n 118 Sn 



exact 


approximate 


exact 


approximate 


2.3571 


2.3571 


2.3311 


2.3311 


3.0115 


3.0116 


2.8941 


2.8942 


3.2495 


3.2496 


3.6361 


3.6361 


4.0229 


4.0229 


4.5070 


4.5070 


4.7795 


4.7807 


4.5569 


4.5577 


5.1974 


5.2010 


5.4049 


5.4069 


5.9277 


5.9320 


5.5611 


5.5620 


6.0754 


6.0790 


5.7648 


5.7673 


6.4924 


6.4968 


6.6544 


6.6556 


6.5964 


6.5984 


6.9403 


6.9429 



Table 2. As 


in table 1, for v 


= 2 states. 




116 Sn 




118 Sn 


exact 


approximate 


exact 


approximate 


2.9832 


2.9839 


2.6957 


2.6958 


3.2691 


3.2694 


3.0218 


3.0218 


3.8398 


3.8403 


4.2904 


4.2905 


4.5816 


4.5817 


4.5047 


4.5088 


5.1152 


5.1263 


4.6708 


4.6720 


5.2917 


5.2973 


5.0742 


5.0742 


5.7828 


5.8034 


5.1998 


5.2109 


5.8861 


5.8907 


5.6664 


5.6704 


6.0401 


6.0456 


5.8045 


5.8047 


6.2093 


6.2439 


6.3618 


6.3651 
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Figure 1. (Colour online) Relative errors in the ground state correlation energy of 
Sn isotopes with 27V active neutrons. Solid lines refer to the procedure of section 3 
at various levels of approximations: the number next to each line indicates how many 
iterative cycles have been carried out in each calculation (see text). The dot-dashed 
line shows the same quantity calculated with the approach of [14] based on independent 
J = pairs . 
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Figure 2. (Colour online) Root mean square values of the relative errors in the 
occupation numbers calculated for Sn isotopes with 2N active neutrons (see text). 
Symbols are the same as in figure 1. 
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Figure 3. (Colour online) Root mean square values of the relative errors in the pair 
transfer matrix elements calculated for Sn isotopes with 2N active neutrons (see text). 
Symbols are the same as in figure 1. 
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Figure 4. Relative errors in the ground state correlation energy of Sn isotopes with 
2N active neutrons. The line labeled with triangles shows the results obtained within 
the approach of section 6 with identical quartets. The line labeled with circles refers to 
the calculations in the PBCS scheme (identical pairs). The line labeled with squares 
shows the best results obtained with the procedure of section 3. 



